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Abstract In the present work we have searched the existence of the late time 
acceleration of the Universe. The matter source that is responsible for the late 
time acceleration of the Universe consists of cosmic fluid with the equation of 
state parameter ui = p/p and uniform magnetic field of energy density pg. The 
study is done here under the framework of spatially homogeneous and anisotropic 
locally rotationally symmetric (LRS) Bianchi-I cosmological model in the presence 
of magnetized dark energy. To get the deterministic model of the Universe, we 
assume that the shear scalar (a) in the model is proportional to expansion scalar 
(6). This condition leads to A = B n , where A and B are metric functions and 
n is a positive constant giving the proportionality condition between shear and 
expansion scalar. It has been found that the isotropic distribution of magnetized 
dark energy leads to the present accelerated expansion of the Universe and the 
derived model is in good agreement with the recent astrophysical observations. 
The physical behavior of the Universe has been discussed in details. 



1 Introduction 

The cosmological observations from type la supernovae [TJ[2], cosmic microwave 
background (CMB) and clusters of galaxies [3] etc, all suggest that the expansion 
of the present Universe is speeding up rather than slowing down. This at once in- 
dicates that the baryon matter component is about 4% of the total energy density 
and about 96% of the energy density in the Universe is invisible which opposes 
the self attraction of matter and causes the observed expansion of the Universe 
to accelerate. This acceleration is characterized by the negative pressure and the 
positive energy density and hence violates the strong energy condition. This viola- 
tion gives a reverse gravitational effect. Due to this effect, the Universe gets a jerk 
and the transition from the earlier deceleration phase to the recent acceleration 
phase takes place [3]. 
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In physical cosmology and astrophysics, the simplest candidate for the dark en- 
ergy (DE) is the cosmological constant (A) [5l[6]. However, it needs to be extremely 
fine-tuned to satisfy the current value of the DE density, which is a serious problem 
7 . Alternatively, to explain the decay of the density, the different forms of dy- 
namically changing DE with an effective equation of state (EoS), u = p/p < —1/3, 
were proposed instead of the constant vacuum energy density. There are several 
other possible forms of DE such as quintessence (cu > —1) [8], phantom (uj < —1) 
[9] etc. It is observed that current cosmological data from SN la (Supernovae 
Legacy Survey, Gold sample of Hubble Space Telescope) [TOlfTT] . CMBR (WMAP, 
BOOMERANG) [12l[T3] and large scale structure (Sloan Digital Sky Survey) [H] 
do not support possibility of u; << —1. However, time-dependent DE characterized 
by ui = — 1 and crossing the phantom divide line even now is a favorable candidate. 
The limit obtained from observational results coming from CMBR anisotropy and 
galaxy clustering is —1.44 < uj < —0.92 with 68% confidence level |14U15j . 

Under the above circumstances, it is observed that in recent years Bianchi uni- 
verses have been gaining an increasing interest and tremendous impetus of obser- 
vational cosmology. In connection to the WMAP data [T5l[T6l[l~7] it is now revealed 
that the standard cosmological model requires a positive and dynamic cosmological 
parameter that resembles the Bianchi morphology 18, 19, 20,21 . According to this, 
the Universe should achieve the following features: (i) a slightly anisotropic spe- 
cial geometry in spite of the inflation, and (ii) a nontrivial isotropization history of 
Universe due to the presence of an anisotropic energy source. The anomalies found 
in the cosmic microwave background (CMB) and large scale structure observations 
stimulated a growing interest in anisotropic cosmological model of Universe. Here 
we confine ourselves to model LRS Bianchi-I whose spatial sections are flat but the 
expansion or contraction rate are direction dependent. For studying the possible 
effects of anisotropy in the early Universe based on the present day observations 
many researchers [22,23,24,25,26,27,28 have investigated Bianchi type-I models 
from different point of view. We notice that, in connection to anisotropic equa- 
tion of state for dark energy, several works are now available in the literature [29, 
30,31 . Some Authors 32,33,34,35,36,37,38 have studied anisotropic DE models 
even with constant deceleration parameter (DP). 

We would like to further mention that unlike the FRW model this Bianchi- 
I type model describes a different kind of Universe in which the scale factor is 
not restricted to be the same in each direction. In the present work, following 
King and Coles [39], we assume a large scale homogeneous magnetic field which 
is responsible for the anisotropy in the flat Universe. It is expected that such a 
magnetic field will impose a single preferred direction in space. Therefore, under 
the influence of this directional magnetic field along the field lines the anisotropy 
of the spacetime will be axisymmetric. However, it is argued by King and Coles 
39j that the observed level of isotropy in the CMB places tight constraints upon 
the strength of any Hubble scale magnetic field. The upper limit on such a field 
is found to be of the order of 10 -9 Gauss [30j. King and Coles |39j also state 
that these and other similar limits quote the adiabatically expanded, present-day 
equivalent very weak value of the field strength which equate to a much stronger 
field at very early times. 

It is a obvious question then - does the above size of the magnetic field sat- 
isfy all observational constraints at the different scales? An extensive field survey 
provides that (i) the lowest measured intergalactic fields and close to the observa- 
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tional upper limits via Faraday rotation measurements [41,42,431, may well be of 
cosmological origin; (ii) a similar protogalactic field strength is inferred from the 
detection of fields of order 10~ 6 Gauss in high redshift galaxies [33] and in damped 
Lyman alpha clouds |45] , (iii) primordial nucleosynthesis constraints only limit the 
equivalent current epoch field to be less than about 1CP 7 Gauss [46 , a value that 
is only slightly stronger than the dynamical constraint at nucleosynthesis [47,48, 

Regarding the impact of magnetic field it is argued that in the early times, the 
magnetic field had the significant role on the dynamics of the Universe depending 
on the direction of the field lines [501139] . Several authors have used Bianchi models 
to investigate the influence of the magnetic field on the evolution of the Universe. 
It is worth noting that there has been some work on magnetic fields in Bianchi I 
models in the past [49,51,52,53 . Milaneschi and Fabbri [51_ studied the anisotropy 
and polarization of CMB radiation where as, Jacobs [39] explored the effect of a 
uniform primordial magnetic field. Both the investigating groups used Bianchi-I 
model of the Universe. Jacobs [39] argued that in the early stages of the evolution 
of the Universe, the magnetic field produced large expansion anisotropics during 
the radiation-dominated phase whereas it has negligible effect during the dust- 
dominated phase. 

In the relatively recent works, it is seen that King and Coles [39| have used 
the magnetized perfect fluid energy-momentum tensor to discuss the effects of 
magnetic field on the evolution of Universe. Sharif and Zubair [54] have studied 
dynamics of Bianchi-I universe with magnetized field of anisotropic dark energy. 

In the present work, however, we present a magnetized dark energy model 
with time varying DP in LRS Bianchi-I spacetime. The investigation is organized 
as follows: The metric and field equations are presented in section 2. Section 3 deals 
with the exact solutions of field equations and physical behavior of the model. The 
comparison between distance modulus (p,) of derived model and observational p(z) 
is presented in section 4 respectively. Finally the results are discussed in section 5. 



2 The Metric and Field Equations 

We consider the LRS Bianchi type I metric of the form 

ds 2 = -dt 2 + A[t) 2 dx 2 + B{t) 2 (dy 2 + dz 2 ) . (1) 

where, A and B are functions of t only. In the limit where A(t) = B(t), the metric 
equation (TJJ) reduces to flat FRW metric. Here the geometry of space-time (1) 
is represented by two equivalent transverse directions y and z and one different 
longitudinal direction x, along which the magnetic field is oriented i. e. A = A x 
only. 

The Einstein's field equations, in the units 8irG = c = 1, read as 

y,i< /;. (2) 

where Tj is the energy momentum tensor cosmic fluid and it is given by 

Tj = diag[p - p B , -up + p B ; , -up - p B , -uip - p B ], (3) 
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where p is the energy density of cosmic fluid, pg is the energy density of mag- 
netic field, u) = p/p is the EoS parameter of cosmic fluid and p is the pressure of 
the cosmic fluid. It is important to note here that the EoS parameter (ui) is not 
necessarily constant [55] . 

It is to note that the energy momentum for cosmic fluid given in Eq. (3) is 
a system of perfect fluid and magnetic field in a comoving coordinates i.e. T- = 
T^(PF) + Ej. Here, Ej is the electromagnetic field. We choose the magnetic field 
along 2>direction. In our model, the electromagnetic field tensor F lJ has only one 
non-vanishing component, viz., F yz = constant. Therefore, for the electromagnetic 
field tensor Ep one gets the following non-trivial components, E t = E% = —Eljj = 
-E z z =p B - 

Therefore, the Einstein's field equations ([2]) for the line-element ([1]) reduce to 
the following system of equations 

2§ + §2 =-^P + PB, (4) 
A B AB 

A + B + AB=- UJp - pB > (5) 

- +2 —= P + PB . (6) 

Here, and in what follows, over-dots indicates differentiation with respect to t. The 
energy conservation equations Tj^ = 0, related to the two equations for cosmic fluid 
and magnetic field |39j , are as follows 

p + 3(l + u)pH = 0, (7) 

pb = Jr- (8) 

Here ft is positive constant and H is the mean Hubble parameter, which for LRS 
Bianchi-I space-time can be defined as 

with a being the average scale factor of LRS Bianchi type-I model and can be 
expressed as 

a = (AB 2 )i. (10) 

It seems that any radiation field with the evolution law expressed in Eq. (8) 
would work. However, one may raise a tricky question that why does this field only 
scale with the scale factor B1 It can be verified that if one writes the conservation 
equation for cosmic fluid and the magnetic field separately, then one can get Eq. 
(8) which contains only scale factor B. 

The spatial volume (V) is given by 

V = a 3 = AB 2 . (11) 

The expansion scalar (6), the shear scalar (a) and the mean anisotropy parameter 
(Am) are defined as 

6 = 3H =J+ 2 §> ( 12 ) 
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where //^(i = 1,2,3) represent the directional Hubble parameters in the direction 
of x, y and z respectively. 



3 Solutions to the Field Equations 

In order to solve the field equations completely, we constrain, the system of equa- 
tions with proportionality relation of shear (<r) and expansion scalar (9). This 
condition leads to the following relation between metric potentials 

A = B n , (15) 

where n is the positive constant. For anisotropic model n ^ 1. 

The motivation behind the assumption, given in the Eq. (|15|) . can be explained 
with reference to the work of Thorne [37] . The observations of the velocity-redshift 
relation for extragalactic sources suggest that Hubble expansion of the Universe is 
isotropic today within approximately 30 percent or to put more precisely, redshift 
studies place the limit a/H < 0.3 on the ratio of shear (cr) to Hubble constant (H) 
in the neighborhood of our galaxy today [561157] . In this connection it is also to 
be mentioned that Collins et al. [58 have pointed out that for LRS type spatially 
homogeneous spacetime, the normal congruence to the homogeneous hypersurfaces 
satisfy the condition a/9 as constant which leads to the assumption A = B" . 

Eqs. flU, © and ([15]) lead to 



(17) 



The general solution of the Eq. (|16[) is given by 

dB 



where £ = 71i ^_ n - ) , k = ni(n — 2), n\ = -^p^ and c\ is the constant of integration. 
Hence the spacetime ([1]) is reduced to 

ds 2 = ^ j-, -+B 2n dx 2 + B 2 (dy 2 +dz 2 ). (18) 



nlB-^ + l)^ B ^ +cl 



After using the suitable transformation of coordinates, B = T, the above model 
(1181) transforms to 



(19) 



6 




0.25 0.5 0.75 



1.25 1.5 1.75 2 2.25 2.5 
T 



Fig. 1 Plot of deceleration parameter q versus time T. 




Fig. 2 Plot of EoS parameter ui versus time T. 



It is to be noted that Eq. (|17[) indicates the explicit dependence of B on t, i.e. 
B = B{t). However, one can not solve Eq. (|17[) in general. So, in order to solve the 
problem completely, we have to choose either B or n in such a manner that (|17p be 
integrable. It can be easily checked that different suitable values of n generate the 
numerical solution of Eq. (|17[) . But we are looking for a physically viable model of 
Universe and this prompt us to consider the transformation B = f(t) = T, which 
gives the time dependent DP. With this transformation Eq. (|17|) leads t = constant 
for T = 0. 

Now, Eq. (H5l) yields 

A = T n . (20) 

It is important to note here that the derived model recovers isotropy with n = 1. 
However, when we put n = 1 in Eq. (|16[) . it leads to a singularity. Therefore one 
can not choose n = 1 to describe the feature of Universe in the present model. 

The physical parameters such as the directional Hubble's parameters (H x , H y , H z 
the average Hubble parameter (H), the expansion scalar (6), the spatial volume 
(V) and the scale factor (a) are, respectively given by 



H x = nniT" 



-(n+2) 



2/ 



(fc + 1) 



(fc+i) 
T "i +ci, 



(21) 
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H V = H Z = mT 



-(n+2) 



21 



(fc + 1) 



T "i +ci, 



(k + iy 



3H = T 



-(M- a )J_M_ T ^ + 



(fc+1)" 



(22) 

(23) 

(24) 
(25) 



The value of DP (5) is found to be 



fc+1 



(26) 



(27) 



The sign of q indicates whether the model inflates or not. A positive sign of 
q corresponds to the standard decelerating model whereas the negative sign of q 
indicates inflation. The recent observations of SN la [1,2, 59_ and CMB anisotropics 
[60 disclose that the expansion of the Universe is accelerating at present and it was 
decelerating in past with a transition redshift about 0.5. It is therefore expected 
a signature flipping in the DP for the Universe which was decelerating in the 
past and is accelerating at the present time [61]. In standard cosmology the DP 
naturally evolves with time just because there are many fluids, and whichever 
takes over determines it at any given time. However, recently Saha and Yadav [62] 
presented an anisotropic DE model with time varying DP. In the present work, 
Fig. 1 depicts the variation of DP versus cosmic time as representative case with 
appropriate choice of constants of integration and other physical parameters. 

The shear scalar (<r) and the mean anisotropy parameter (A m ) are given by 



18T 2(n+2) 



( 3nni ) ( 3ni ) 



2( 



1 



3nni 



+ 2 1 



(fc + 1 



3ni 



(fc+i) 
■T -1 +ci 



(28) 



(29) 



The energy density of the cosmic fluid (p), the EoS parameter (cj) and the energy 
density of magnetic fluid (pb) are found to be 



^271+2 



2n_ fr 

rp rp^ 



2n 1 £T~' 1 + (3T - 2)T 2n+1 - /3T 2 



,2 f^L T ^ 



+ Cl 



1 + 3*- 



PB — yj- 



(30) 

(31) 
(32) 
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Fig. 3 Plot of shear scalar a versus time T. 



The parameters H, 6 and a start off with extremely large values and continue to 
decrease with the expansion of the Universe whereas the spatial volume (V) grows 
with the cosmic time. Fig. 2 shows the variation of EoS parameter (cj) versus 
cosmic time for accelerating phase of the Universe as a representative case with 
appropriate choice of constants of integration and other physical parameters. 

In the derived model the shear is limited to about 10 -4 (Fig. 3) which is in 
fair agreement with the work of Adamek et al [63] . 



4 Some Observational Constraint 

In this section we follow the maximum likelihood approach under which one min- 
imizes \ 2 an d hence measures the deviations of the theoretical predictions from 
the observations. 

Let us now provide the scale factor a and redshift z are connected through the 
relation 

a 



(1 + *)' 



(33) 



where ao is the present value of scale factor. 

Combining Eqs. (|23[) . (|26|) and (|33|) . one can easily obtain the expression for 
the Hubble's parameter (H) in terms of redshift parameter (2) as follows 

3(n+4) 

H = H Q {1 + z)-^^. (34) 

Here Hq is the present value of the Hubble's parameter. Note that this relation is 
obtained by omitting the constant of integration. 

It is necessary for the investigation of type la supernovae to explore DE and 
constraint the models. Since SN la behave as excellent standard candles, they can 
be used to directly measure the expansion rate of the Universe upto high redshift, 
comparing with the present rate. The SN la data gives us the distance modulus 
(fi) to each supernova as 

H = 5log w D L (z) + fio , (35) 
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0.12 



Fig. 4 The distance modulus fi as function of redshift (z) to the derived model compared 
with SN la data from Amanullah et al. I64l . The observational fJ.(z) data points are shown 
with error bars and the solid dots corresponds to distance modulus of derived model. 



where Dl = H °^ L is the Hubble-free luminosity distance and fio is the zero point 
offset, defined as 



f^o = 5log 10 



Mpc 



+ 25. 



Inserting Eq. (|36[) into Eq. (|35[) . we obtain 



[i = 5log w 



( d L 
\Mpc 



+ 25. 



The luminosity distance (d^ is calculated by 

djj = ri(l + z)a . 



(36) 



(37) 



(38) 



For the determination of n, we assume that a photon emitted by a source with 
co-ordinate (r,t) and received at a time to by an observer located at r = 0. Then 
we determine n from following relation 



n 



■° dt 
a 



(39) 



By solving the Eqs. (|37p ~ (|39p and (|26p . one can easily obtain the expression for 
distance modulus (/i) in the term of red shift parameter (z) as 



fj, = 5log 10 



3H (n+2)k (i + z y 



21 <«+2)<= 
k+1 



(1 + 2)— -1 



+ 25, 



(40) 
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Table 1 Comparison of the results of present model with the observational data 





Redshift(z) 


Supernovae Ia(fi data ) 


Our model (/J. th ) 


0.0087 


30.35lg;» 


32.07943 


0.0104 




32.46627 


0.0135 


33 in+ ' 22 

3>5.1U_ 22 


33.031435 


0.0172 


33.93lg_ 


33.555816 


0.0245 


34 33 +0 A 


34.320871 


0.0285 


34 61 + - 23 


34.64755 


0.0306 




34.8010 


0.0365 


35.181°;" 


35.18139 


0.0453 


35.9llg;i? 


35.64668 


0.0499 


35.84_g;iS 


35.85475 


0.0529 


36.09l°;« 


35.98026 


0.0589 


36.35 q y 


36.21104 


0.0603 


36.33io.l7 


36.26146 


0.0627 


36.75l°;« 


36.34521 


0.0701 


36.71±g;_| 


36.58437 


0.0746 


36.86_°;1° 


36.71760 


0.0786 


36.76_°;i«j 


36.82936 


0.0876 


07 tr+0.19 

°''°-0.19 


37.06104 


0.1009 


01 ■ lo -o.n 


37.36251 



where Hq is in the unit of Km s Mpc . 

In the present analysis, we use 19 data set out of recently released 38 data 
set of SN la in the range 0.0015 < z < 0.12, as reported by Amanullah et al. |64| 
(Table 1). In this case Xsn nas been computed according to the following relation 



X 2 SN=A-^+log 1Q (j^ =97.772, (41) 



where 



38 ( data _ ,,th\ 2 
\L (1 j 



E 



38 ( data _ th\ 
yH 1 ) 



M 1 



The comparison between distance modulus fi of the derived model and the 
observational fi(z) SN la data as reported by Amanullah et al. [6_] can be seen in 
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Fig. 4. The observational fj,(z) data points are shown with error bars and the solid 
dots corresponds to distance modulus of the derived model. It is observed that the 
derived model is fit well with SN la observation (see Fig. 4). Note that the best fit 
value of distance modulus is fi(z = 0.0365) = 35.18139 with Xsn = 97.772 and 
the reduced \ 2 value is x^jvAdegree of freedom)=2.64. 

5 Results and Discussions 

In this paper, we have investigated magnetized DE model under the assumption 
that a oc 9 in Bianchi-I spacetime. Under some specific choice of the parameters the 
present consideration yields the time dependent DP and EoS parameters. It is to be 
noted that our procedure of solving the field equations are altogether different from 
what Sharif and Zubair (5j have adopted. The derived model starts expanding 
with Big Bang singularity at T = and this singularity is point type because 
the directional scale factors A and B both simultaneously vanish at T = 0. The 
dynamics of DP parameter yields two different phases of the Universe. Initially DP 
is evolving with positive sign that yields the decelerating phase of the Universe 
whereas in the later times it is evolving with negative sign which describes the 
present phase of the acceleration of the Universe. Thus the derived model has 
transition of the Universe from the early deceleration phase to the later acceleration 
phase which is in good agreement with the recent observations [J. 

The distance modulus of the derived model fit well with the observational 
n(z) values (see Fig. 4 and Table 1) which in turn imply that the derived model 
is physically realistic. It is important to note here that in the absence of the 
magnetic field only the anisotropic distribution of DE leads the present acceleration 
of the Universe [37] while in the presence of magnetic field along with the isotropic 
distribution of DE describes the dynamics of the Universe from Big Bang to the 
present epoch. Thus the magnetic field isotropizes the distribution of DE which 
signifies the role of magnetic field. Hence from the theoretical perspective, the 
present model can be a viable model to explain the late time acceleration of the 
Universe. In other words, the solution presented here can be one of the potential 
candidates to describe the observed Universe. 

It is worth noticing that we have presented solutions of a magnetized Bianchi-I 
Universe like King and Coles [39] . However, their work is concerned with vacuum 
energy. In our case, we have studied the Universe consists of cosmic fluid with 
equation of state w = p/p. Therefore, our study is more general than King and 
Coles [32]. 

In the present work we observe an interesting feature that the magnetic pres- 
ence isotropises the expansion in the model. However, it is not clear, in particular, 
what property of the magnetic field is responsible for the isotropisation. Is it the 
additional energy density of the field, is it the magnetic pressure, or maybe the 
tension? It can be observed that the solution of Eq. (16) does not exist for n = 1. 
Therefore, in the derived model, one can not choose n = 1. To obtain an explicit 
solution of Eq. (16), one can choose /3 = which leads to pg = and the isotropic 
distribution of cosmic fluid. This can also be observed from the LHS's of Eqs. (4) 
and (5), which are same for n = 1. As a result, one should get pg = 0. This means, 
there is no contribution of magnetic field and the model becomes an isotropic 
FRW model. The well known criteria for isotropization are A = ^Sjfr — 1 — > 
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and E 2 = \AH 2 — 1—^0, where A, E 2 are average anisotropy and shear respec- 
tively, and H 2 is the directional Hubble parameter. Obviously, these criteria are 
valid for large physical time. In the present model (Eqs. (21) - (24)), one can easily 

find out that for large time, A — )■ 2< ^ l +2 ]' 2 an< ^ ^ ~ * This immediately implies 
that for n — >• 1 our model turns out to be an isotropic FRW model. Therefore, 
n = 1 is the condition of isotropy in the absence of magnetic field and the presence 
of magnetic field has constraint on n = 1. This seems to contribute the magnetic 
field which resembles with the initial anisotropy of the Universe. With the passage 
of time magnetic field decreases and becomes negligible at late time to approach 
towards isotropy. 

Our solution, in the present investigation, shows that the constrained equation 
of state of p evolves in time, but this is not all that is needed when comparing 
with observations. This explicitly means that it is not exactly known about the 
generic fluid p, and how does it describe and include the known cosmological 
history where an early radiation domination gave way to matter domination. In 
particular, fitting the H{z) evolution law derived from SN observations is only one 
of the many pieces of information which need to be used. A thorough discussion 
of all basic cosmological constraints is beyond the scope of this analysis in the 
present paper and is awaiting for a future project. 
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